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ABSTRACT 


This thesis is comprised of four independent computer 
programs and their related operating instructions. Each of 
these progrags focuses on a particular facet of tactical 
missile design. The topics covered include guidance and 
trajectory calculations, rocket motor propulsion sizing, 
warhead design, and aerodynamic coefficient determination. 
The programs are developed from accepted mathematical proce- 
dures and ara tailored to optimize operator interaction for 
educational purposes. This thesis is intended to be uti- 
lized as a supplement to the thesis Tactical Missile Concep- 
tual Design by D.&R.Redmon, Naval Postgraduate School, Sep- 
tember 1980. 
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I. Dg 9) 


The programs contained in this thesis were assenbled 
expressly to supplement the work of Dan Redmon in his the- 
sis, Jactical wissile Conceptual Design. Iwo of the 
programs, LPATH, the trajectory model, and LAERO1, the aero- 
dynamic coefficients program, originated in Redmon's thesis 
and were expanded/modified for use on the Naval Postgraduate 
School's new IBM 370 computer system. The other two pro- 
gtams were generated for this thesis and utilize the proce- 
dures and principles outlined by Redmon. 

The specific intention of these programs is to provide 
students of tactical missile engineering and design with a 
method of solving complex mathematical routines rapidly and 
interactively. Each of the programs request data which are 
likely to be used as design parameters for the topics con~ 
cerned. The programs also allow repeated operation with 
input alteration capability, allowing the user to personally 
optimize his design. This approach was chosen to allow stu- 
dents to understand the relationships various input parame- 
ters have with the final solutions. 


II. RAJ ORY MODELS 


A. DESCRIPTION AND ORIGIN OF PROGRAM 

This program applies the principles of missile guidance 
laws to the terminal phase (the last 5 to 10 seconds) of a 
missile trajectory in order to determine the maximum normal 
acceleration on the missile for a given scenario. Of the 


three general guidance law categories, pursuit guidance is 


not included in the program capability. It has been found 
that pursuit guidance invariably produces a tail-chase situ- 
ation, greatly reducing an anti-air aissile's effectiveness 
against maneuvering targets of similar speed characteris- 
tics. Line-of-sight guidance and proportional navigation 
guidance are both options of the prograa. 

Figure (II-1) shows a typical encounter geometry as 
required for this program. The encounter is considered to 
cccur entirely within a two- dimensional plane. No differ- 
entiation is required or assumed concerning the orientation 
of the encounter plane. The plane may be at any angle to 
the horizontal as desired by the program user. The refer- 
ence direction is an arbitrary choice by the program user. 
The angles shown are positive in value, however the program 
dees not require positive angles. If TAL were 150 degrees, 
it could aiso te entered as -210 degrees. The IRA term 
represents the initial range to the target. 

Tangential velocities {air speeds) of the missile and 
the target are considered by the program to be constant 
throughout the problea. Since the program concerns itself 
with only tha final moments of a trajectory, this is a rea- 
sonable consideration. Target normal accelerations, when 
specified by the user, are also held constant throughout the 
prcblem for the same reason. The missile normal 


Program Variables 

LOS’ (line of Sight angle) 
target alpha) 
missile alpha) 
target speed) 
misSile speed) 
initial range to. 
target from missile 
launch _ point) 
(initial range to 
Missile from 
launch point) 


Angles are shown as 
positive quantities 


Reference direction 


Pigure (II-1). Encounter geometry 
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( acceleration, however, is a function of the guidance laws 
! and will vary according to the resulting flight path 


4 requirements. A constraint on the encounter is that the 
| missile normal acceleration must be zero at the start of the 
1 preblem. 


The program "flies" both the missile and the target as 
Sisple points in space with no consideration given to aero- 
dynamics. The missile will always strike the target dead 
ceoter when it is given the proper speed advantage for the 
encounter since there is no provision for statistical miss 
“Jt analysis. The unbounded nature of the missile normal accel- 
§ eration allows the missile to turn as rapidly as necessary 
a to hit the target. 
This program analyzes the given encounter by time incre- 

ment calculations. as is the case with all integrations 

conducted by incremental steps, the accuracy of the results 

is a function of the increment size. The results will 

become increasingly accurate as the tigze increment is nade 

swaller. However, as the time increment decreases in size, 
= the length of the output becomes increasingly longer. The 

user must balance the desire for accuracy against the amount 
| of time he wishes to spend on the computer terminal. } 
: The primary output is a tabular listing of the aissile 
| and target coordinates at each time increment. The coordi- 
a nate frame is established within the encounter plane with 
{ the abscissa oriented along the reference direction. The 
| problem stops once the missile has passed its closest point 
of approach to the target. Output then provides the time to 
intercept from time of problem initiation and the maxiaua 
} acceleration the missile was required to endure. A Versatec 
k. plot of overlaying successive encounters is an optional 


output. 
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This program originated as two seperate BASIC programs 
woitten by Redmon (Ref. 1} for use on the HP 9830 desktop 
calculator. It was subsequently translated into FORTRAN IV 
for use on the IBM 370 computer systen. 


B. USER INSTRUCTIONS 

If it is desired to abort the operation of this prograg 
prematurely, two methods are available. When the progras is 
waiting for data entry, press [ENTER]. when the program is 
not waiting for data entry but is processing, type "HX" and 
press (ENTER]. Both actions will return the terminal to CMS 
mode. 

When the screen becomes full, or "MORE...." appears in 
the status area, clear the screen by pressing and 
Simultaneously. At several points in the prograa, 
the user will be directed to "CLEAR SCREEN AND INPUT 0". 
This is to provide proper positioning of the output on the 
screen for ease of reading. If any other symbol than 
a/b should appear in the lower left of the screen, ; 
press and continue. 

1. Turn the terminal on with the red Switch. 

2. When the large "NPS" appears after about 30 seconds, 
press [RESET], then press [ENTER]. 

3. When "CP READ" appears in the status area on the 
lower right of the screen, type "L nnnnP", where onnn is 
your 4-digit user number, then press (ENTER]. 

4. You will now be asked for your password. Type it in 
(the characters will not appear on the screen), then press 
(2uTER]. 

5. Your personal file must contain a PROFILE EXEC rout- 
ine with the appropriate Fortran GLOBAL statement. If it 
does not, type "GLOBAL TXTLIB FORTMOD2 MOD2EEH", then press 


(EnTERy . 
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6. Type "LINK TO xxxxP 191 AS 192 BR", where xxxx is 
the 4-digit user number for the project file, then press 
[Eurea]. 

7. You will now be asked for the project file password. 
Type it in (the characters will not appear on the screen), 


then press [{ENTER|. 
8. Type “ACCESS 192 BY" and press |ENTER]. 


9. Press [KIT] and simultaneously to clear 
screen. 
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pees ios 7 aa a " 
Ae han. 2m ne es A 
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ae 10. Type "LPATH" and press (ENTER). 
| 11. Input the following data as it is requested on the 
el screen by typing the data and then pressing [ENTER}]. Ensure 
i that the data is input as either decimal or integer as spe- 
= cified. The terms in parenthesis below are the program 
‘ variable names. 
adl 
= Parameter Units Value range 
Trajectory option (TITLE) none 0=Line-of-sight 
\ 1=Proportional 
= navigation 
a 2=Both 
Time increment (DEL) seconds Larger than .0001 
: times the problea 
time 
Navigation constant (NAV) none 2.5 to 6.5 
LOS Angle (LOS) degrees 0.0 to 360.0 
| Initial target range meters Larger than the 
{ froma the launch site (IRT) misSile range 
¢ Target speed (TSP) n/sec Larger than 0.0 
Target noraal n/sec/sec positive is to 
{ acc@leration (TAC) target's left 
pt 
Initial target angle to degrees 0.0 to 360.0 the 
| line of sight (TAL) 
Be Missile speed (MSP) a/s ec Larger than 0.0 
Initial aissile raoge meters Such that impact 
frca@ launch site (IM) occurs in less 


Pe than 10 seconds 
12. If you desire to rerun the problea, or want to run a 

new problea, answer the questions appropriately when asked 

by the terminal after the execution of the current problea. 
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13. To veceive the printout and plot of your encounters, 
answer "no" to rerunning or restarting the problem when 
asked by the terminal and follow the directions presented on 
the screen. 

14, Upon completion of the program, type "LOGOFF" and 


press | ENTER]. 
switch. 


15. Turn the terainal off with the red 
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Figure (II-2). Non-maneuvering crossing target 


Rt=10000 meters 
Rm=9000 meters 
60=30.0 degrees 
Otoz=-90.0 degrees 
at=0 

Vt=200 meters/second 
Va=800 meters/second 
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Table (II-1) is the corresponding computer output for 
this encounter. As indicated, the aissile maaxinum normal 


acceleration is 
am= ~32.00 m/sec/sec or -3.26 g's. 
Figure (II-3) is the Versatec plot of the engagement. 
2. Examole II-B. ; j ve 
crossing Target 
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Figure (II-4%). Maneuvering crossing target 


Rt=10000 meters 
i Ra=9000 meters 
0230.0 degrees 

‘ Ato=-90.0 degrees 

| at2=156.8 m/sec/sec (16 g's) 
Vt=200 meters/second 
Va=800 meters/second 
Navigation constant=3.06 
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Table (II-~2) is the computer output for this encounter. 
The missile maximum norrgal acceleration is 
am= 283.94 m/sec/sec or 28.95 g*s. 
Figure (II-5) is the Versatec plot of the engagement. 
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D. PROCEDURAL FLOWCHART 


bee ets Input data ; 


Establish initial angular 
relationships and detérmine 


_anitia eof 
X-Y coordin tes for missile 
and target 


SERIE iiaianaaiaa 


Increment iteration counter 


Have you exceeded 9999 itera- 
tions? 


Sani: Saee 


|_____erint ercor_sessage | 


Print error message 


w-n- PERE SESE MEH. 


[ew 


Calculate angular relation- 
ships between missile and tar- 


get 


Calculate noraal acceleration 
on aissile 
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e2 ee wae2ee ae wen bwe eee 2 aee2ew eee woe 


Determine range between ais- 
Sile nad target 


Seenahiiaias ania 


Calculate new X-Y coordinates 
for missile and target based 
on speeds. 
and angular relationships 


Increment time | 


Is range greater than zero? 


eee SL Mens Pr 
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PO ox Print output | 
Do you wane _iust line- of- 
: ght? 
§ yes 
END i 
"Establish initial Sree 
relationships and detérmine 
; initial Xx-y 
coordinates #t aissile & tar- 
"4 ge 


Increment iteration counter | 


~Tagzensh’_tenptesuneee._ 


\ Gave you exceaied 9999 dtera= 


<u you exceeded 9999 itera- 
ons 


Ent siea 


Print error message 
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E. PROGRAM CHANGES 


1. Language Translation. 


The two BASIC pro cams contained in Referense 1 were 

translated into standard FORTRAN IV. 
2. program Condensatiog. 

The two individua programs were combined to form a 
single integrated routine which allows the user *o choose 
either of the two guidance laws or both for a given encoun- 
ter. The two original programs were meshed such that only 
the actual guidance law algorithms are separate routines, 
all input and output rout nes are now common. 

3. Input and Output Facility. 

The data input instructions were modified to maxin- 
ize user facility on the I3M 370 System 3278 terzinals. 
Coapletion of data input is now followed by data feedback 
for user verification prior to actual program execution. 

Data output has been expanded to provide data deliv- 
ery to three destinations; to the user's terminal for obser- 
vation, to a printer file for hardcopy duplication 9f the 


terminal display, and to a plot data file for subsequent use 
by the plot program. The iestinations are options chosen by 


the user for each execution 9f the program. Up to nine 
different problems can be printed andi plotted each time the 
program is entered. 
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The program can now be rerun multiple times without 
exiting and re-entering each time. The user has a choice of 


either rerunning the same problem or initiating a new prob- 


lem completely. 
4. Plot Program. 

PATHPLOT FORTRAN was developed to produce a Versatec 
plot picture of the encounter. It will produce a single 
plot each time the program is entered and will plot up to 9 
engagements in an overlaying manner. This format was 
adopted to allcw comparisons of successive input data nodi- 
fications. 

5. Data Qvercapacity Check. 

If the user initiates a problema requiring more than 
9999 iterations, the program will stop. The user will be 
notified of the error and given the opportunity to rerun the 
problen. 

6. {pitial Missile acceleration. 

Initial missile acceleration was removed as a user 
input variable and established as zero. Due to the mathe- 
matical nature (i,e., no physical constraints) of the pro- 
gram, any “wrong” accelerations of the missile in the ini- 
tial state were immediately "corrected" by the algoritha. 
The model is better served by providing no initial accelera- 
tions. 

7. Theta Angles. 

Both the target and the missile theta angles (the 
angles between the velocity vectors and the reference direc- 
tion) were removed as user input variables. The program now 
calculates the theta angles from other input data, reducing 
redundancy and possible contradiction of input data. 

3. Initial conditions Peg spective. 

Originally, the gwissile guidance calculations 

started at t=0. Specifically, the anchor point for the line 
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of sight solution was the point where the missile commenced 
the problem, whether or not that was on the launch site. 
This produced a situation removed from reality where the LOS 
anchor point should be at the fire control location, usually 
at or near the launch site. A similar, though less pro- 
nounced, condition existed for the proportional navigation 
sclution. The program was modified to provide proper posi- 
tioning of the external guidance reference. As a result of 
the modification, an additional output is the correct lead 
angle for tha missile at the start of the problem. This 
angle is based on the assumption that neither the target nor 
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the missile have maneuvered prior to time t=0. 
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IIIT. WARHEAD DESIGN 


A. DESCRIPTION AND ORIGIN 
j This program develops a warhead using the same nethods 
| as presented by D Redmon (Ref. 1]. dowever, its capability 
| is somewhat greater and it applies the relationships in a 
Slightly different manner. This program was written for 
this thesis in FORTRAN IV for use on the IBM 370 computer 
system 3278 terminals. It is configured expressly for oper- 
| ator-computer interaction. 
It starts with an initial input of data from which a 
table of- fragment initial velocities is generated and pre- 
Lo sented to the user. Prom this, the user selects values for 
fragment mass and impact velocity and another table is gen- 
erated and displayed. From this one, the desired probabil- 
ity of hit given a detonation is selected and the final 
solution is produced. At various points during the opera- 
, | tion of the program, the user will have opportunities to 
alter or revise selected parameters. 
The progras is limited to a cylindrical warhead with 
either a solid or hollow core. The fragments: are required 
to be square in shape and are sized by the program. Figure | 
| (III~1) illustrates the shape of the warhead and the loca- ) 
tion of various input and output quantities. \ 
Initially, the target altitude is used to determine the 
atmospheric density, temperature and speed of sound. These 
values are, in turn, used to calculate required initial vel- | 


ecities for the fragments. The program is preloaded with 
various values for the fragment mass and impact velocity, 
he which are used to generate the table of initial velocities. 

These velocities come from the following relationship: 


Z ¥ (hit) =Vifexp (-ks) J 
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where V(hit) is the impact velocity, Vi is the initial vel- 
ocity, s is the kill radius, m is the fragment mass, Pa is 
the atmospheric density, A is the plan area of the fragment, 
and Cd is the drag coefficient of the fragment. 

The ballistic limit velocity is calculated for the vari- 
ous presized fragments and is provided as a reference when 
choosing an appropriate impact velocity. The ballistic 
liait velocity is that velocity at which one aalf of the 
fragments will penetrate the target's skin and the other 
half will not. The empirical relationship, developed by A. 
E. Fuhs (Ref. 6], presents the ballistic liait velocity as a 
function of the fragment size to skin thickness ratio. dis 
function dealt with steel fragments impacting an aluminua 
plate. His results were modified to qualitatively account 
for different skin materials and fragment densities. 

Next computed is the fragment spray angle and the criti- 
cal aiss distance. The spray angle is a function of the 
initial velocity, the detonation velocity and the warhead 
length. The critical awiss distance is defined as the range 
where the fragment spray exactly covers the entire target. 
The critical miss distance is used by the program to sepa- 
rate the two functions which determine the average number of 
hits received by the target. The program assumes the target 
is always centered within the fragment spray and aligned 
perpendicular to the central ray of the spray. 

A selection of warheads is then produced, one for each 
of a greloaded set of Pd's (probability of a hit given a 
detonation) to provide the user with a Pd versus warhead 
weight/size trade off comparison. This sizing process is 
based entirely on the following relationship: 


Pd=1-exp (-a) 
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shewn by Redmon (Ref. 13, is a function of the cube of the 
warhead radius. The user then chooses a desired Pd which, 
in turn, produces the final warhead sizing. 


- where a is the average number of hits. fhe variable a, as 
| 


{ Values that are shown in tabular form for user selection 
and input into the program are not limited to those dis- 
played. Any value in between the displayed values or con- 
pletely outside of the value range may be used. The one 
exception is that Pd can never be selected to be larger than 
-999 and may even, if forced by the program limitations, be 


4 ee 


required to be lower. Since ultimately in this algoritha, 
Pd is a function of the warhead radius, the maxiaunm value 
for Pd may te reduced in order to keep the radius within the 
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original input parameters. The user is notified if this 


at 


condition occurs. 
i Useful reference information for some common explosives 
and case materials is contained in the following tables. 


5 Table (III-1). Characteristics of common explosives 


Osiv bepsaly paced) ) v t/sec) 

RDX - 0596 9300. 26837. 

HMX - 0665 10230. 29934. 

: FETN - 0625 9300. 27231. 

: TETRYL - 0585 8200. 25755. 

| COMP B - 0607 8800. 25722. 

OCTOL - 0650 9500. 28356. 
1 


Table (III-2). Densities of common case saterials 
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CTH 


FRLEN 


Individual 
fragment 


RO=sDW 


RC=4sDW-CTH 


| 
| 
| 
| 
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Core 


Explosive 
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Pigure (III-1). Form of warhead as used by the program 
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B. USER INSTRUCTIONS 

If it is desired to abort the operation of this program 
prematurely, two methods are available. When the program is 
waiting for data entry, press [ENTER]. When the program is 
not waiting for data entry but is processing, type “HX" and 
| press (ENTER]. Both actions will return the terminal to CMS 
mode. 

when the screen becomes full, or "MORE...." appears in 
the status area, clear the screen by pressing and 
simultaneously. At several points in the prograa, 
the user will he directed to "CLEAR SCREEN AND INPUT 0". 
This is to provide proper positioning of the output on the 


aL 
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screen for ease of reading. If any other symbol taan 
" eo should appear in the lower left of the screen, 


press |RESET] and continue. 
switch. 


a 


1. Turn the terminal on with the red 
2. When the large "NPS" appears after about 30 seconds, 
press [RESET], then press [ENTER]. 
t | 3. when “CP READ" appears in the status area on the 
a lower cight of the screen, type "L nnanP", where anno is 
your 4-digit user number, then press [ENTER]. 

4. You will now be asked for your password. Type it in 
{the characters will not appear on the screen), then press 
(estes). 

5. Your personal file must contain a PROFILE EXEC rout- 
ine with the appropriate Fortran GLOBAL statement. If it 
does not, type "GLOBAL TXTLIB FPORTMOD2 MOD2EEH", then press 
[ERRER} . 

6. Type "LINK TO xxxxP 191 AS 192 RR", where xxxx is 
the 4-digit user number for the project file, then press 
(emrER] . 

7. You will now be asked for the project file password. 
Type it in (the characters will not appear on the screen), 


then press |ENTER|. 


ee . 
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8. Type “ACCESS 192 BY and press | ENTER|. 
9. Press and {CLEAR| simultaneously to clear 
1 screen. 


| 10. Type "LBOMB" and press [ENTER]. 
11. Input the following data as it is reguested on the 

screen by typing the data and then pressing [BNTER]. Ensure 

that the data is input as either decimal or integer as spe~ 

cified. The terms in parenthesis are the program variable 

- nages. Input the following variables as decimal values: 


. Explosive density (XDEN) lb/cu.in 
+ Explosive Gurney constant (GC) ft/sec 


4 Explosive detonation velocity (VD) ft/sec 
- Case material density (CDEN) lb/cu. ft 
Desired kill radius (RKILL) feet 
7 Warhead diameter (DW) inches 
Warhead Llength-to-diameter ratio (LTD) 7 
Target length (LT) feet 
Target mean width (WT) feet 
| Target vulnerability, P(k/h) (PKH) 
Target altitude (ALT) feet 
Target skin thickness (TTH) inches 
Target skin material (MAT) 1.0=Aluainun : 
| 2.0=Piberglass/Kevlar 
{ 3.0=Steel , 


12. After entering the: above data, you will be presented i 
with an initial velocity table built around your desired 

| kill radius. The initial velocities will be listed as a 

function of fragment sass and impact velocity. Also pro- 

: vided will be the ballistic limit velocities for each of the 

fragaent aasses. Input the following parameters as decimal 


values: 


grains 


Fragment mass (IFMLB(1)) 
Impact velocity (VHIT(1)) ft/sec 


ts 2 nop ee 


13. You will now be presented with a shopping list of 


warheads developed for a range of Pd's. The warheads are 
described by the following quantities; 


Warhead weight in pounds (total weight) 
Case thickness in inches 

Core diameter in inches 

Number of fragments from the warhead 
Number of fragments hitting the target 
Edge length of the fragments in inches 


Input the following parameter as a decimal value: 
Desired probability of a hit given a detonation (PDF) 


14. If you desire to rerun the problem, or want to run a 
new problem, answer the guestions appropriately when asked 
by the terminal after the execution of tke current problen. 

15. To receive the printout of your encounters, answer 
"no" to rerunning or restarting the problem when asked by 
the terminal and follow the directions presented on the 
screen. 

16. Upon completion of the program, type "LOGOFF" and 
press [ENTER]. 

17. Turn the terminal off with the red switch. 


C. EXAMPLE PROBLEMS 

1. Example Zii-a 

It was desired to build a warhead which would kill a 

typical cruise missile flying at 100 feet altitude. The 
warhead was required to have 4 lethal radius of 50 feet with 
a Pd of 0.98. TNT was selected for the explosive and steel 
was chosen for the case material. The diameter of the mis-~ 
sile was 13.5 inches. 
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Table (III-3) outlines the input parameters. fTatkle 
(III-4) shows the program output. After the initial veloc- 
ity table was displayed at the terminal, the fragment mass 
was chosen to be 130 grains and the impact velocity was 
selected to be 2500 feet per second. After the Pd table was 
displayed, 0.98 was entered as the kill probability. 
2. Example [iI-8 


f A warhead was required which would kill a typical 
| 
| 


ny Yano 


manned aircraft at 30000 feet. A lethal radius of 75 feet 
was specified. The warhead was limited in weight to 50 
{ pounds and needed to have a core diameter of at least 5 
e inches. The diameter of the missile was 10.0 inches. HMX 
, was chosen as the explosive and depleted uranium as the case 
material. 
Table (III-5) outlines the input parameters. Table 
(III-6) shows the program output. After the initial veloc- 
ity table was displayed at the terminal, the fragment mass 
was chosen to ke 210 grains and the impact velocity was 
selected to be 5000 feet per second. After the Pd table was 
displayed, it was determined that the fragments were too 
large and the fragment mass was then reduced to 100 grains. 
The impact velocity was also reduced to 3000 feet per sec- 
: ond. When the Pd table was redisplayed, 0.995 was chosen as 
) | the desired kill probability. 
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INPUT DATA FOR EXAMPLE ITI-A 


TABLE (ITI-3). 
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INPUT DATA FOR EXAMPLE III-B 


TABLE (III-5). 
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OUTPUT DATA FOR EXAMPLE III-B 


TABLE (III-6). 
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TABLE (III-6). 
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D. PROCEDURAL FLOWCHART 


Calculate characteristics of 
atmosphere at altitude 
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Calculate initial velocities 
for all combinations of preset 
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velocities and fragment 
weights at the kill radius 
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Do you want t9 change the kili 
radius? 
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Input new kill radius 


Input desir i fragment mass 
and impact velocity 


“Calculate the critical frag- 
ment spray radius 
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Calculate charge-to-mass ratio 
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Determine warhead solution for 
specified pd 
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‘4 IV. PROPULSION MOTOR SIZING 


Aa. DESCRIPTION AND ORIGIN OF PROGRAM 

This program provides a method for the preliminary siz- 
ing of a solid propellant rocket motor for a boost-sustain 
trajectory of a tactical missile. The analytical method was 
developed by Redmon (Ref. 1] and was expanded with the addi- 
tion of material from Platzek [Ref. 2] and Hill [Ref. 3]. 
4 The program was written for this thesis in FORTRAN IV for 
sf use on the IBM 370 computer system. Essentially, the analy- 
! sis consists of sizing the booster motor from differential 


Sg, 


velocity and acceleration requirements with limitations 
imposed by the physical dimensions of the missile. The 
bocster is at all times considered to be a core-burning 
motor. The sustainer motor calculations are controlled ty 
the maximum cange specified by the user and by the solution 
Z| of the booster motor. The sustainer can be either a core- 
burning or an end-burning motor. 
The rocket motor configuration is assumed to be either 
an integral bocster-sustainer motor as shown in Figure 
(IV¥-1) ox a staged motor as in Pigure (IV-2). The booster 
: | and the sustainer always burn exclusively, or, in other 
| words, one is not burning while the other one is. The 
nozzle half angle is specified by the user consistant with 
Space available in the missile. If a staged motor is uti- 
lized, both nozzles will have the same half angle. 

The booster calculations start by determining the amount 
of thrust needed to boost the total weight of the missile to 
its cruise velocity at the prescribed acceleration. From 
this, the necessary amount of propellant is estimated and 
the process is iterated to account for the decreasing mass 
situation. The chamber pressure is initially estimated by 
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Pigure (IV-1). 
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Figure (IvV-2). 
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{ miniaizing the motor weight to propellant specific impulse 
ratio, as presented by Redmon [Ref. 1]. The next step of 
the program is to size the nozzle using the now known ini- 
tial pressure ratio. Once the ideal nozzle is developed, 
the chamber pressure is raised or lowered as necessary in an 
attempt to cause the nozzle exit diameter to exactly match 
the booster diameter. However, if the integral notor option 
is used, the program will drive the nozzle as small as it 
can without violating one of the following limits in order 
to increase the probability that the sustainer will operate 


eR aS Shek ee 
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properly. The iterative process has two limits: an abso- 
lute smaximaum of 2000 PSI chamber pressure [Ref. 2] and a 
DBinimum of 1000 PSI if the pressure had previously been 
higher. The solution can be less than 1000 PSI if the pres- 
sure remained below that level throughout the problem. Also 
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limited is the exit pressure to ambient pressure ratio. At 
the low end, it is limited to 0.4 to prevent flow separation 
| in the nozzle. At the high end, it is limited to 2.5 to 
4 prevent excessive underexpansion and loss of physical real- 
ity in the program results [Ref. 4). The burn rate is 
initialized at a potential maximum of 1.25 inches/second 
(Ref. 2] and is allowed to decrease to arrive at a compati- 
ble burn area and web thickness combination. 
: | The sustainer motor, in the integral motor case, is vir- 
F tually a continuation of the booster solution. The initial 
\ thrust requirement is determined by increasing the cruise 
| speed drag to account for speed loss through maneuvers. [t 
| 


is also then refined for weight change if any climbing or 
4 diving is reguired to reach the target altitude. The nozzle 
is the same one as developed for the booster except that the 


* throat area will be expanded as a result of the erosion 
effect. From the new area ratio, a pressure ratio is deter- 


os mined. The chamber pressure and thrust coefficient are then | 
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iterated until a steady chamber pressure evolves to provide 
the required thrust. If at any tige it drops below 125 FSI, 
the program will stop since this is considered a mnininua 
chamber pressure for proper propellant combustion. The exit 
pressure to ambient pressure ratio remains subject to the 
same restrictions. The burn rate starts at 0.45 inches/sec- 
ond {Ref. 2} and is decreased to provide an acceptable web 
thickness and burn area. The solution can be either an end 
burner or a core burner, depending on the burn area 
required. 

The sustainer for the staged motor is solved in essen- 
tially the same manner as the booster. The two major excep- 
tions are that it can produce an end burner if the burn area 
is saall enough and the thrust required is based on the 
cruise drag and not the velocity to be gained. The chamber 
pressure is limited to an absolute maximum of 800 PSI anda 
minimum of 250 PSI if the pressure had previously been 
higher (Ref. 2]. The Same pressure ratio restrictions apply 
and the burn rate starts at 0.45 inches/second. 

Other motor-nozzle combinations can be created from the 
results of this program on a first-order approximation 
basis. Figures (Iv¥-3) and (IV-4) illustrate two possible 
methods for separating the motor nozzles without resorting 
to staging. Although these two methods will probably have 
the same weight disadvantage that staging does, they can be 
packaged in to a smaller volume of space. A hint for "crea- 
tive construction" is to rerun the problem after sizing the 
motors using a smaller missile diameter to force the nozzles 


to a smaller size. 
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B. USER INSTRUCTIONS 

If it is desired to abort the operation of this program 
prematurely, two methods are available. When the program is 
waiting for data entry, press (enreR]. When the program is 
not waiting for data entry but is processing, type "Hx" and 
press [ENTeR]. Both actions will return the terminal to CMS 
ucde. 

When the screen becomes full, or “MORE...." appears in 
the status area, clear the screen by pressing and 
Simultaneously. At several points in the program, 
the user will be directed to “CLEAR SCREEN AND INPUT 0". 
This is to provide proper positioning of the output on the 
screen for ease of reading. If any other syabol than 
« AR should appear in the lower left of the screen, 
press [RESET] and continue. 

1. Turn the terminal:on with the red switch. 

2. When the large "NPS" appears after about 30 seconds, 
press [RESET], then press [ENTER]. 

3. When "CP READ" appears in the status area on the 
lower right of the screen, type "L nonnP", where nnna is 
your 4-digit user number, then press [ENTER]. 

4. You will now be asked for your password. Type it in 
(the characters will not appear on the screen), then press 
; 

5. Your personal file aust contain a PROFILE EXEC rout- 
ine with the appropriate Portran GLOBAL statement. If it 
does not, type "GLOBAL TXTLIB FORTMOD2 MOD2EEH", then press 
[Eres]. 

6. Type "LINK TO xxxxP 191 AS 192 RR", where xxxx is 
the 4-digit user number for the project file, then press 
(esrea]. 

7. You will now be asked for the project file password. 
Type it in (the characters will not appear on the screen), 


then press |ENTER]. 


49 


1 
‘ 
4 
( 


fae 


Rr pect aoe ges a a hk tc rt te a 


8. Type “ACCESS 192 B" and press | ENTER]. 
9. Press and [CLEAR] simultaneously to clear 
screen. 


10. Type "LPROP" and press [ENTER]. 


11. Input the following data as it is requested on the 
screen by typing the data and then pressing [enTeR]. Ensure 
that the data is input as either decimal or integer as spe- 
cified. The terms in parenthesis are the program variable 


Dages. 
The following variables are required inputs for both 


actor option problems. 


Motor option (IMOTOR) Q=integral motors 
1=staged motors 

Launch altitude (LALT) feet 

Launch weight (WL) pounds 

Launch valocity (¥ BI) feet/second 

Launch elevation angle (ELB) degrees 

Boost acceleration (A) gravities 

Cruise velccity {V BF) feet/second 

Cruise velocity drag (DRAGS) pounds 

Maxigum cange to target (R) nautical miles 

Maxigum target altitude (IALT) feet 

poceiticriabises (I SPB) seconds 

Booster propellant : 

density (DENSB) pounds/cu. inch 

Boosters exhaust. 

specific heat ratio {G B) 


Sustainer prope. tant 
specific ispulse 


Sustainer propellant 
density i! 


(I SPS) seconds 


(DENSS) pounds/cu. inch 


Sustainer exhaust 
specific heat ratio (GS) 


Nozzle half angle (ALN) degrees 


Bs es ee a ee a a ae 


The following variables are required inputs for the 


integral motors option only. 


Nozzle design altitude (ALTBN) feet 


Nozzle erosion rate (EB) inches/second 
Missile diameter (DB) inches 

Case yield strength (YIELD) PSI 

Case density (DENSC) pounds/cu. inch 


. ta 
The following variables are required inputs for the 


staged motors option only. 


Booster design altitude (ALTBN) feet 


Booster diameter (DB) inches 

Booster case 

yield strength (YIELD) PSI 

Booster case density (DENSC) pounds/cu. inch 
Sustainer design 

altitude (ALTS N) feet 


Sustainer diameter (D) inches 


Sustainer case 
yield strength (YIELD) PSI 


Sustainer case density ‘=(DENSCS) pounds/cu. inch 


12. This program will cue the user when the input paraa- 
eters have dictated a scenario which either cannot be 
achieved within reality or produce less than optinun — 
requirements on the propulsion system of the missile. They 
are not definitive and are only intended to make the user 
aware of a situation which may need correction. The follow- 
ing is a list of available cue messages with short 
definitions. 

“SUSTAINER NOT CALCULATED SINCE THE BOOSTER BURNOUT 
RANGE EXCEEDS THE DESIGN RANGE." This can result frea 
entering an extremely short range for the missile, or it 
can be Caused by a very low accéleration requirement. 
Paebeeeene Geek Gan Obeecet THE post ee NOZZLE DESIGN 
PREVENTS SUSTAINER OPE ATT ON. MMEND STAGING OR 
THDEDENDENT? HOZZLES." athis Sccurs “OnLy when us sng the 
° 


integral motors option. Then scenario described the 
program can cause the booster nozzle to be too large to 
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maintain the proper chamber pressures when the motor has 
shifted to sustainer operation. This usually occurs 
when a large acceleration is demanded but the thrust 
required fcr the cruise trajectory is small. 


“BOOSTER BURN RATE WAS ADJUSTED TO CORRESPOND WITH THE 
WEE THICKNESS." This occurs guste oftea and Siepi) 
indicates that the burn rate was decrease rom it 
petential physical maximua of 1.25 inches/second. 


“THE BOOSTER NOZZLE DESIGN WAS RESTRICTED DUE TO THE 
MISSILE DIAMETER,” and "THE SUSTAINER NOZZLE DESIGN WAS 
RESTRICTED DUE T6 THE MISSILE DIAMETER." The nozzle was 
not able to be designed for optimum pressure ratios at 
the mid point of the boost trajectory. Usually, the 
exit diameter is solved larger than the missile diameter 
and is subsequently reduced to fit. 


“BOOSTER NOZZLE DESIGN [S NOT OPTIMUM DUE TO EXCESSIVE 
BOOSTER CHAMBER PRESSURES," and "SUSTAINER NOZZLE DESIGN 
IS NOT OPTIMUM DUE TO EXCESSIVE SUSTAINER CHAMBER 
PRESSURES." If the nozzle cannot be downsized without 
exceeding Prosser thresholds (2000 PSI for the booster 
and 800 PSI for the sugtainer) the chamber pressure is 
held just kelow the pressure freshold and the nozzle 
area ratio will be adjusted to allow the nozzle to fit 
in the missile. 


“THE SUSTAINER MOTOR HAS A CORE-BURNING GRAIN." The 

required burn area for the sustainer was too large tco 

permet an end burning grain with a properly realistic 
urn rate. 


"“SUSTAINER BURN RATE WAS ADJUSTED TO CORRESPOND WITH THE 
WEB THICKNESS." This indicates the sustainer burn rate 
was lowered from a potential maximum of 0.45 inches/sec- 
ond to provide a proper web thickness. 


"THE SUSTAINER MOTOR HAS AN END BURNING GRAIN." The 
required burn area for the sustainer was small enough to 
petmwit an end burning grain. The burn rate is then 
adjusted to correspond with the nonreduceable burn area. 


“REESTIMATION OF LAUNCH WEIGHT IS REQUIRED FOR Tyee. 
MISSILE PERFORMANCE PARAMETERS." The scenario described 
to the program produced a motor whose weight is either 
larger than 75% of the total or less than 25% of the 
total. 
"ENLARGEMENT OF DIAMETER IS RECOMMENDED DUE TO A VERY 
HIGH LENGTH-TO~DIAMETER RATIO FOR THE MOTOR.” This cue 
indicates the length-to-diameter ratio is greater than 
- .Other components of the missile will make the mis- 
silefs overall length~to-diameter ratio even larger. 
13. Immediately after completion of the solution, the 
ercgram will ask if you want to receive a hardcopy printout 
of that particular solution. A “yes" answer stores that 
scluticn in a file for retrieval by the user when he fin- 


ishes with the prograa. 
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14. If you desire to rerun the problem, or want to run a 
new problem, answer the questions appropriately when asked 
by the terminal after the execution of the current problen. 

15. To receive the printout of your encounters, answer 
“no" to rerunning or restarting the problem when asked by 
the terminal and follow the directions presented on the 


screen. 
16. Upon coapletion of the program, type "LOGOFF" and 


press {ENTER(. 
17. Turn the terminal off with the red switch. 


C. EXAMPLE PROBLEMS 


1. Exaaple [y-a. t motors, common nozzle 


The following parameters are input for the integral 


motor example: 


Launch altitude 35.0 feet 

Launch weight 1000.0 pounds 
Launch velocity 0.0 feet/second 
Launch angle 60.0 degrees 
Average acceleration 30.0 g*s 

Cruise velocity 4000.0 feet/second 


Drag at cruise velocity 1500.0 pounds 
Maximum range 20.0 miles 
Pinal (target) altitude 50000 feet 
Booster propellant 


specific impulse 260.0 seconds 
Bocster propellant . 
density 0.075 lbs/cu.inch 


Booster exhaust specific 
heat ratio 1.244 


Sustainer propeitant 
Specific impulse 


Sustainer propellant 
ensity 


210.0 seconds 


0.065 lbs/cu.inch 


Sustainer exhaust 


specific heat ratio 1.270 
Nozzle half angle 20.0 degrees 
Nozzle design altitude 0.0 feet 
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OUTPUT OF EXAMPLE IV-A 


TABLE (IV-1L). 
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7 Nozzle erosion rate 0.001 inches/second 
7 Missile diameter 10.0 inches 
| Yield strength of 
¥ case material 180000.0 PSI 
Density of case : 
material 0.2662 lbs/cu.inch 
The solution for this problem is presented in Taktle 


Ans 


(IV¥-1). 
2. Example Iv-8. staged motors, separate nozzles 


The following are input for the staged motor prblem: 


4 Launch altitude 35.0 feet 

i Launch weight 2200.0 pounds 

# Launch velocity 0.0 feet/second 

5 Launch angle 30.0 degrees 

. Average acceleration 25.0 g's 
Cruise velocity 2200.0 feet/second 
Drag at cruise velocity 1000.0 pounds 
Maximum range 50.0 niles 


| Pinal (target) altitude 75000 feet 
| Booster propellant 


specific impulse 250.0 seconds 
Booster propellant : 
density -0647 lbs/cu.inch 
: Booster exhaust. 
| specific heat ratio 1.225 
Sustainer, propellant 
a specific impulse 205.0 seconds 
Ei Sustainer propellant ' 
density 20625 lbs/cu.inch 
; Sustainer exhaust 
specific heat ratio 1.257 
a Nozzle half angle 15.0 degrees 
Booster design altitude 0.0 feet 
, Booster diameter 14.5 inches 
Yield strength of : 
booster case 180000.0 PSI ; 
“ Density of booster . 
case material 0.2662 lbs/cu.inch 


Sustainer design altitude 0.0 feet 
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GUTPUT OF EXAMPLE IV-B 


TABLE (IV-2). 
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Sustainer diameter 14.5 inches 


Yield strength of 
sustainer case 180000.0 PSI 


Density of sustainer : 
case material 0.2662 lbs/cu.inch 


output for this example is provided in fable 


a: dey 
gM  ee 


' 
| 
1 
: 
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D. PROCEDURAL FLOWCHART 


Input data i 


SS SEeniid GAGE 


esewvwoeowoew we emotes war eee 2 ewe owe 


Calculate average drag of 
boost phase based on cruise 
speed drag and initial and 
final speeds 
“Estimate thrust required for 


boost and calculat? booster 
burn time 


Saaaniiid iiaaeeaeen 


Determine average altitude 
during boost and calculate 
ambient pressure 


Refine thrust required for 


boost and determine yoight of 
booster propellan 


we2ee ce oe nee wee ewe aww oon 


|"~“betermine a we wae awe amerwoow | 


Determine anass flow rate 


~~-2esensdne. sft ee Ente 


ew oer woe aewoeew ewe Swe eae er wooe aawewoere 


Calculate pressure ratio based 
on optimum, hamber pressure 
to case weight relationship 


Saisie Viaabneihaan 


|______ButHs2000 pst | 


PLIM=2000 PSI 


Calculate thrust coefficient, 
throat dianetar and exit 
ameter 


Is exit diameter smaller than 
base diameter? 


2 e2ee @e eae eaee ge a 2eae2eeens wea wo 


NO 


| _inetease chasber pressure __| 


Increase chamber pressure 


Lower Chamber pressur?2 
Increase exit pressure 


Is (exit pressure/ambient 
pressure) greater than 0.4? 


Peace epee eee et 


NO 


awweew ee woeee eee bow wee @ eee eee a2oee 


Is (exit pressure/ambient 
pressure) less than 2.5? 


Set exit diameter equal to 
base diameter 


Calculate area ratio 


Oe wee Oe eee eee @ ee eee es eee ew eae 


Enererer nr, Spann Saeuias 


Is chamber pressure less than 
PLIM? 


J xo 


nosso agate seasaate 


Lower chaaber pressure 


Sei sewcece seek cacceccemeoe a VS 


Is this an integral motor? 


NO 
Is the nozzl2 exit diameter 
4 larger than 95% of the base 
dianeter? 
Wo 


Is the chamber Rrosegre less 
h PSI? 
~ “Lower PLIM 
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erase eee eee 


Set burn rat equal to 1.25 
and calculat2 rasulting web 
thixskness 


required web thickness 


Are the two web thickness 
values equal? 


Determine burn area and | 


i. 


aii as il 
ee 


Calculate di aensions of the 
booster and nozzle 


ke oat a bs ite 


wwe wee oe we www bow owe ow ewe oe ew oe Yes 


< Is this a staged motor? 


wo 
Calculate required sustainer 
burn time 


Determine average altitude 
during cruise and calculate 
ambient pressure 


Calculate thrust required to 
offset craise drag and 
Calculate weight of sustainer 
propellant 


Determine mass flow rate and 
nozzle area ratio 


[Calculate wach of the exhaust 


sense eee ESS ee tee eeneees 


Calculate nozzle pressure 
ratio and estimated thrust 
coefficient 


Calculate sustainer chamber 
pressure 


SEcnnnil Sica 


<s chamber pbessure less than 


wo 


| ee ee @ ee eee owe bee wee eo owes eee oo 


Refine thrust coefficient 


Does thrust coefficient equal 
estimated thrust coefficient? 


eceweee were we 2 om we we Oe@e eee eee ae 


Set burn rats equal to 0.45 


Set PRED Eat SEUsE 2 02. 


Determine required burn area 


“sasatatay WEEE aT Gee | 


Is the burn area less than the 
motor cross sectional area? 


Is the web thickness less tha 
86% of the motor radius? 


NO 


( 
| 
| over the burn rate | 


Lower the burn rate 


eee 2e2e ee eee “weaece 222 2 e@e eee eo 


Set burn ar a equal to motor 
cross-section area and lower 
burn rate to sorresponad 


~—-#aen ESE Co ceetseeet 


Establish dinensions of sus- 
tainer aotor 
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SCE Sle vel, SS ae 


para 


Oo Se 


+: 


‘“- 


oy salt” < Salieri. ne 


Calculate sustainer burn time | 


Determine average altitude 
during cruise and calculate 
ambient pressure 


ee were wee eee eo oe 


Calculate thrust required to 
offset cruise drag and 
calculate weight of sustainer 
propellant 


| Determine aass flow rate 


~-teteen ine. 15s flee tase. 


Calculate pressure ratio based 
on optimum hamber pressure 
to case weight relationship 


Siasianeaniaid Geieasiainiion 


PLIM=800 PSI 


[eeeeennaa PEED G0 Bets | 


Calculate thrast coefficient, 
throat diasaeter and exit 
diameter 


SEARLS shh, Sia ene eee 


s exit diameter smaller than 
base diameter? 


eo wwe ew eee ceoae oe 


Seaee © ee 2e 2e @ we ga oe ete wee es eeen eee 


Ie Sse Bw wwBenaehaes ears ear aeen oe 


i aa | 
na 
a 
fa 
r) 
| 
o 
+O 
n 
pp 
a0 
n 
w 
a 
an 
i=] 
wh 
@ 
=) 
ry) 
nH 
wa 
oO 
4H 


Lower chasber pressure 
Increase 2xit pressure 


tee ee 


“Is (exit pressure/ambient_ Yes 
pressure) greater than 0.4? 
Tuo 


Is (exit pressure/ambiant 
pressurej Lass than 2.5? 


Ses Dees ewe Pos ener eeweeees oe 


Set 2xit diameter equal to 
base diameter 


Calculate? area ratio 


Is chamber pressure less than 
PLIM? 


cae: f" ee 


Lower chaaber pressure 


Is nozzle exit diameter Larger 
than 95% of the base diameter? 


Tne 
Is the chamber pressure less 
than 250 PSI? 
wo 


ene ewe we we aew 5. eee | 


Low2r PLIM 


22 2 ee ee ee ween we eww wen wom ew eee ew on 


| Determine case thickness 


Calculate burn area 


arisiaiiis Wisma 


Is burn area smaller than 
motor cross-section area? 


Calculate required web thick~ 
ness and web thickness 


{ 
| 
{ 
| 


ase I TRE RR IE Li A ANT he ciate hd a ha Png thtce tints abi ts a cide 6. m 


el ene lem am. 
» 
“4 
0 
ct 
or 
Dp 
~~ 
o 
Q 
(= 
p 
~ 
bod 


ae een 
Lower burn rate 


Set burn ar@a, equal to motor 
cross-section area and 


1 
| calculate corresponding burn 
| rate 
; { 
aan: en ee eee wo ew wwe howe waco eee ewe won 
| Calculate dimensions of sus- 
+ tainer motor and nozzle 
| nn Mice auseoees (3) 
| rr Ark ree ee ee ae Pee Ser re 
a Set all sustainer values to 
Z zer 


Print output | 


waa EEE 


ee eee 


eras 
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A. DESCRIPTION AND ORIGIN 

This program is the current edition of a program which 
originated at the Naval Ship Research and Development Center 
in 1971 (Ref. 5}. The program was written in FORTRAN II for 
use on the IBM 7090 digital computer. It was developed asa 
method for predicting the static, longitudenal aerodynamic 
characteristics of typical missile configurations with the 
ccntrcl surfaces in a "plus" attitude. The original program 
computed the aerodynamic characteristics for missiles oper- 
ating at angles of attack up to 180 degrees. The effects of 
control surface deflections for all modes of aerodynanic 
ccotrcl are taken into account. The method was based on the 
then well known linear, nonlinear crossflow and slender tody 
theories with empirical modificaticns to provide the high 
angle of attack capability. 

The program was modified and, presumably, updated in 
19574 by F. A. Kuster, Jr., of the Naval Air Development Cen- 
ter. In 1980, the program was modified for use on the Naval 
Postgraduate School IBM 360 computer system by D. Redmon 
CRef. 1J. The current version of the program was modified 
fcr use on the new Naval Postgraduate School IBM 370 compu- 
ter systema. It has been expanded to provide graphical pre- 
sentation of the output data. 

It must be eaphasized at this point that the current 
program edition is not in a completely finished state. 
Scmewhere in the history of the program after its initial 
estapiishment on the IBM 360, errors were introduced during 
the modifications. At present, these errors do not prevent 
the use of the program and the output data is considered to 
be correct for trend-observance purposes. 
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Specifically, the prograa does not produce any drag-rise 
phenogwena for either the wings or the tails when Cd is 
cbserved as a function of Mach number. Additionally, the 
decline of the drag coefficient above Mach 1.0 is not smooth 
or aS prolonged as is found experimentally. It is very fFro- 
bable that these two failings of the program are linked to a 
ccmmon error inserted accidentally in the process of tailor- 
ing the program for use on the IBM 370. In order to 
temporarily smcoth over the graphical discontinuities, expo- 
nential decay functions were inserted. They are clearly 
marked in the the prograg listing for removal when the pro- 
gram is corrected. 

The input to the program is composed of a detailed list- 
ing of the dimensions of the missile to be analyzed. The 
current version of the program will consider a missile which 
has four symmetrical wings and four symmetrical tails. The 
missile may be either canard or tail equipped and either 
wing or canard or tail controlled. The program assumes the 
centrol surface is the tail, however, the input data is 
“aislabelled" to produce the proper configuration. For ins- 
tance, if the gwissile is a wing control aissile, the wing 
data is input as the tail and the tail data as the wing. 

For a canard controlled missile, the canard data is input as 
the tail. Figures (¥-1) and (V-3) Show two typical aissile 


configurations and where the input parameters for the pro- 


gram are obtained. 


Be. USER INSTRUCTIONS 

If it is desired to abort the operation of this program 
prematurely, two methods are available. When the program is 
waiting for data entry, press [ENTER]. dhen the prograns is 
not waiting for data entry but is processing, type "HX" and 
press [ENTER]. Both actions will return the terminal to CMS 
aode. 


wins Chita, ot lg iaalemuleaaapt acai: Bg Nn ane sare gst ts al ne ae Vihear ine ales nit, 


When the screen becomes full, or "MORE...." appears in 
the status area, clear the screen by pressing and 
simultaneously. At several points in the prograan, 
the user will be directed to "CLEAR SCREEN AND INPUT 0". 
This is to provide proper positioning of the output on the 


eae ei, en SiMe se gene 


screen for ease of reading. If any other symbol than 
"GAR « should appear in the lower left of the screen, 


press [RESET] and continue. : 
switch. 


1. Turn the terminal on with the red 


aes 


2. When the large "NPS" appears after about 30 seconds, 


press [RESET], then press [ENTER]. 


3. When "CP READ" appears in the status area on the 


‘qoee8 
‘plain — svete! cane Gen 


Ne 


lower right of the screen, type "L nnonnP", where nnnn is 
your 4-digit user number, then press [ENTER]. 
4. You will now be asked for your password. Type it in 
(the characters will not appear on the screen), then press 
(exrza] 
5. Your personal file must contain a PROFILE EXEC rout- 
oa ine with the appropriate Fortran GLOBAL statement. If it 
does not, type "GLOBAL TXTLIB FORTMOD2 MOD2EEH", then press 
(Earea]. 
6. Type “LINK TO xxxxP 191 AS 192 RE", where xxxx is 
: the 4-digit user number for the project file, then press 
| LETT 
| 7. You will now be asked for the project file password. } 
Type it in (the characters will not appear on the screen), 


rad 


then press [ENTER]. 

8. Type “ACCESS 192 B" and press [ENTER]. 

9. Press and simultaneously to clear 
screen. 

10. Type "LAERO1" and press [ENTER]. 

11. Input the following data as it is requested on the | 
be screen by typing the data and then pressing [ENTER]. Ensure 


* 
(eit. Mag: tees oe 
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that the data is input as either decimal or integer as spe- 


a te 5 Aaa ir ORCS A CE Re EN 


cified. The terms in parenthesis below are the progran. 
variable names. 


Input the following as integer values unless otherwise L: 
noted. The integers aust be two digit integers (three=03). ! 


Cenatrol constant (ICSC) 01-Tail control 
Q2-Wing control 
03-Canard control 


Nose shape (INOSE) Qi-Ellipsoid 
02-Ogive 
03-Cone 


Nuaber of control deflections {F0T) Less Shap 11. : 
You wiil now be asked for e control deflections in 
degrees, decimal values. 


Number of Mach numbers (7M) Less than 25 

You will now be asked for the Mach numbers, decimal 
i values. Each Mach number will produce a separate 
21 table and plot of output data. 


Nuaber of angles of attack (LAL) Less than 25. 
You will now be asked for the angles of attack in 
degrees, decimal values. 


Number of confi urations (HBODT No restrictions 
Bach configaration will restart the program. Only 
‘ the last coOnfiguration will produce the written 


Leek SMe Pe oe 


oF output. 
Wing planform (ISWPW) Qi-Not delta , 
02-Delta 
Wing position (IAFBW) at a after wing. 
01-No body after wing 
Wing sweep constant (ISWEP®) 00-Delta planform or 
pu suee leading 
: edge 
| 01-Swept leading edge 
a Nuaber of wings (NWING) 04 
gi Tail planform (ISWPT) Q1-Not delta 
@ Q02-Delta 
7 Tail position (TAFBT) 00-No body after tail 
Q1-Body after tail 
Tail sweep constant (ISWEPT) 00-Delta planform or 


Unswept leading 
edge ; ; 
Q1-Swépt leading edge 
Number of tails (NTAIL) 04. 
Input the following values as decimal numbers: 
* Wing tip-to-chord ratio (xLAMW) 
Wing leading edge sweep (CLAMW) Degrees 
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ey 


es 


hoa cba ee ae See ee eee 


na 


wing span including body (bw) 
Wing root chord (CROOTH) 

Exposed wing area, 2 panels (S#) 
ding wean geometric chord (XMACW) 


Distance from nose to 
wing leading edge (XWING) 


Wing thickness-to-chord ratio (TOVCW) 


Tail tip-to-chord ratio (XLAMNT) 
Tail leading edge sweep (CLAST) 
Tail span including body (5T) 
Tail root chord (CROOTT) 

Exposed tail area, 2 panels (ST) 


Tail mean geometric chord (XMACT) 


Distance from nose to 
tail leading edge (XTAIL) 


Tail thickness-to-chord ratio (TOVCT) 


Altitude (HT) 

Body diameter (D) 

Missile length (XL) 

Nose length (XLNOSE) 

Distance from nose to CG (XCG&) 
Reference area (AREA) 
Reference length (XREF) 
Engine code (ENGINE) 


Inlet code (ENLET) 


Boat taii angle (BETA) 

Base diameter (DBASE) 
Nozzle exit diameter (DJET) 
Boat tail length (XLABOD) 
Protuberance drag (CDPROT) 


If comparing with experimental 
values, Reynolds number (REFT) 
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Peet 
Feet 
Square feet 
Feet 


Feet 


Degrees 
Feet 

Feet 
Square feet 
Feet 


Feet 


Feet 
Feet 
Feet 
Feet 
Feet 
Square feet 
Feet 


0.0-Turbofan 
1.0-Rocket 


0.Q-Flush 
1.0-Extended 


Degrees 

Feet 

Feet 

Feet 

(Coefficient value) 


(Dimensionless) 


va we 


oe each ae 


ea vd 
ate 


aaa: 


12. If you desire to rerun the problem, or want to run a 
new problem, answer the questions appropriately when asked 
by the terminal after the execution of the current problea. 

13. To receive the printout and plot of your encounters, 
answer "no" to rerunning or restarting the problem when 
asked by the terminal and follow the directions presented on 
the screen. 


14. Upon completion of the program, type "LOGOFF" and 


press [ENTER]. 
15. Turn the terminal off with the red switch. 


C. EXAMPLE PROBLEMS 
Table (V¥-5) identifies the output variables as they 
appear in the output tables. 
1. Example v-4. Jail control missile 
Pigure (V-1) illustrates the missile used in this 
example, The dimensions for this missile and cther input 
parameters are contained in Table (¥-1). The output is 
shown in Takle (V¥-2) and Pigure (¥-2). 
2. Example V-IL. Canard control missile 
Figure (V¥-3) illustrates the conard configuration 
Bissile used in this example. The input data is contained 
in Takle (V¥-3). The output is displayed in Table (V-4) and 
Pigure (V-4). 
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CROOTT 


Pigure (V¥-1). Tail control missile as used in Exaaple V-A 


73 


: 
: 


THE FCLiCHING TABLE CONTAINS THE INPUT DATA FOR 


TAIL CCNTAOL MISSILE 


EXAMPLE V-A. 


CDOONMNOTCONDLOQOTHOOOOrS oo 
BDOOFKFOMDOOOANANEOOONOND ow 
CODTOPFOCOMNOR FF YOMINC FON NAHMSOQOFNOO 
wh TNL oral mat om Al mt af ee teererevstereteeveaet*e teeters veneee 
~ OR ANN AMOOU MOOONO SS DOMOD4D00000 
oO 


io] 
oO 
toa] 


9 
4 
1 
l 


3=C ANARD: 
3=CONES 
LEADING EDGE: 
LEADING EDGE: 
JNCTION} = 


NOSE: 


be) Jeb b- OR 

UUNTegani zz 
mn 

SAOWLD tee MWD PLN ee LNO 
ed DDD ADL NQ eR TQ RDO. ZN DE 
we ASCP qNmYAZIMNAD YE WZ Ay we 
ZW IDA LAUD WOE OIIT A Ne UIA . 
ANA SWSZ FROZE Rr OL DIWYSIY DrwWoOSHM ee LOD Pe et 
eh JD — SON Die MY Diy Ds vei Tt ee NU Foe DN > 
Ae 2LOIOTA £EDIOREMALZTY LOX ALTO WIXS Leas en ze z 
ZYUIIA A & UM (Dette Ih ZI. Li~ wr S A9I0D Tru 
ORS TIO ZF COTOKOYIZ oO X9Z0 oe Few 
Ut TF QWOERRAWMOAEZBEOOR DADRA Omit wwe wwe IIT TW 
Te & QOQAWYN BOM Iw egewWOT EM erator FWY OUW se RO 
S VIL IO AED SOI CPN RI eer ss UT TP ety 61D 
uU U we OdDuwZ jDwe au a2rauy AW WOW Wee 
MUY CoD Irina y 2) JOR OM YG gO ket eee uw 
PAW WS SOE FLA. ZNGW ZN SR NWN DS ID BZ IRANI 
ZUOMwWOIVW ££ DWH Wa Syme NT ee RL a QD 
UZxXDASHAO <0 Zeeyo OFO RRR OF DRS YY MW MUDD 
uw DED ~ zZz~5. ere HP Naneeem 2D ee 
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VZOZ2ADAM BNAY Hd EY £F EVI eM Et OR Mm tee ew IYO 
mt att ae Zee tt KD TI KO ICICI Me TOI KM MU DOIN KU 


ee et ee a et re a er mr cee et er ee tee ee Se ee Se ee Nt et eet See Set PN a Net SP ey See Ne See 


OCATION FROM 


mo aw Oe er = rene RS OR EE, ER Te Ey yy OR Oe Meee eA 
SNAP AVM OM OANA TUNOM DHOOM SIVOE DAOAN DIN OE DFROANMNS 
eh ed med ok ed INIA NII NII NASAL CY AOTC SSPE 


Input data for Example V-Aa 


Table (¥-1). 
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aed Soe 


pei 


CM 


xCP 


XCPA 


CA XCPW XCPT 


CN 


CNWP CNIP CLIO COoID 


col 


CLB 
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o9 sosseosoeses 
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Table (V-2). 
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XCPB 
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Output data plot for Exampl2 V-B 
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Table (¥-5). Output variables of LAERO1 


; AL Angle of attack 
; CLTOT Total coefficient of lift 
= CDTOT Total coefficient of drag 
| CLWP Wing panel soefficient of lift 
CLBW Additional Lift on body due to wing 
i CLTP Tail panel coefficient of lift 
| CL3T Additional Lift on body due to tail 
CLB Body alone Lift coefficient 
¥ CDI Induced drag coefficient 
& CNWP Wing panel normal force coefficient 
‘ CNTP Tail panel normal force coefficient 
3 CLTD Lift coefficient due to tail deflection 
- CDTD Drag coefficient due to tail deflection 
CN Total normal force coefficient 
CA Total axial force coefficient 
XCPW Wing center of pressure 
XCPT Tail center of pressure 
XCP Total missile center of pressure 
CM Total pitching moment about C.G. 
CDOWBT zero lift drag coefficient of wing-body- 
tail combination 
: CDMISC Miscellaneous zero lift drag coefficient 
' | cpDoT zero lift drag coefficient of tail 
a cpow Zero lift drag coefficiant of wing 
| CDOB zero lift decay coefficient of body alone 
CDPROT Drag coefficient of body protrusions 
CDINL Drag coefficient of engine inlet 


ae CDAFT Drag coefficient of boattail region 
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+ VI. CONCLUSION 


e There are many topics which may be the subjects of fol- 
| low on work contained within this thesis. Although the four 
i programs have been installed on the IBM 370 computer systen, 
these four alone do not fully satisfy the original goal of 
this work: Provide a computer supplement to the Tactical 
Missile Conceptual Design handbook. Numerous additional 
focal algorithas are utilized in the design handbook which 
deserve the attention of a programmer. Of immediate inter- 
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est are the areas concerning radar or infrared guidance 

systems, baseline configuration modelling and weight distri- 

bution, and initial control and lifting surface design. 

Each of these topics san be programmed to provide missile 

design students interactive learning tools when coupled with 

the design handbook. 

The most urgent follow on work to this thesis is the 

7 restoration of the program LAERO1 toa reliable, useful pro-~ 
gram. The program was modified and set up on the IBM 370 
computer system during the period immediately following the 
system's installation at the Naval Postgraduate School. As 

| could be expected, the computer suffered many and varied 

1 grcwing pains in its early life. As a result of this, or of 

the human manipulation expanding the capability of the pro- 

gram, the effectiveness of LAERO1 was substantially reduced. 

Work involving the other three programs would involve 

simply expanding their capabilities. The trajectory models 

program, LPATH, presently considers only two guidance laws: 

line-of-sight guidance and proportional navigation guidance. 

Other guidance laws which can easily be included in the pro- 

gram include pursuit, beam rider, and combinations of diffe- 

rent laws. It might also prove useful to be able to 
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( simulate the entire missile trajectory but still only output 
the terminal phase of the encounter. Another option would be 
te provide the target with a controlled trajectory instead 
cf the constant acceleration condition now imposed. 

The propulsion sizing program, LPROP, should have the 
various nozzle options incorporated into the program so that 
it isn't necessary to manually juggle the program output. 
Other booster-sustainer grain configurations could be 
explored, such as the booster grain being cast within the 
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ccre of the sustainer, or even a motor with only a single 
gtain. Another suggestion for the convenience of the pro- 
glam users is to institute a shopping list of availlable 
preopellants and their characteristics into the progran. 
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A similar list of available explosives and case materi- 
als could be put into the warhead sizing program, LBOMB. 


we. 


These shopping lists would provide ready access to reference 

information and, at the same time, decrease the number of 

data values to be manually input into the computer. Since 

| the current program is limited to cylindrical warheads, an 

| area of expansion would be the flexibility of warhead 

styles, such as curved, shaped charge, continuous rod, etc. 

At present, the programs are somewhat hindered by the 

mechanics involved in producing the printouts and the plots. 

Due to the results of tailoring a program to be interactive, 

often it is required to completely exit a program before 

output can ke received. Subsequent design iterations | 
i 
| 
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require re-entry into the program, which produces a certain 

justifyable annoyance to the user. Additionally, the compu- 

ter center has instituted a new policy of cancelling any 

jobs with duplicate job names, which can be severely irri- 

tating and cumbersome to the persons running the plot rout- 

ines contained within this thesis. However, the computer 

ot center has developed procedures which have the potential to 
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alleviate both of these problems. According to Volume 13, 
Number 4 of the Computer Center Newsletter, CMS commands can 
now be invoked from within a FORTRAN program. The print and 
plot operations presently contained within executive rout- 
ines can now be placed directly within the source programas. 
This will remove all current restrictions placed om the nua- 
bers of printouts received per session and will label each 
plot with the user’s job name and not the project's job 
nage. 
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APPENDIX A. TRAJECTORY MODELS PROGRAM PLOWCHART 
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MAIN PROGRAM 


nesta awe waoore oe 
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[ TBACK=0 ] 


ee Sea! ee 


Input data: TPITLE,DEL,NAV,LOS 
IRT, TSP, TAC, TAL, MSE, IRA 


[call screen 
“EPRO=1_ 
TLOS=1 


manana pee 


Playback input data to terminal 
for verification 


<Is data correct? 


MTA=LOS+MAL 
: LOS=DETORA(LOS 
TALSDETORA (TAL 
TTH=DETORA(TTH) 
| MAL©OETORA (MAL | 
MTH=DETORA (NTH 
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T3=TAS/TSP 
P3=TSP*sSin(A1) /R1 


A3=T3-P3 
R3=TSP*cos(Al 
R4=85 P*cos(A2 


D6=R4*TSP*sin (A1) 
Donner teecnecto Tiny 
D8=—-R3*MSP*sin (A2 

DO=R1*MSP*cos (A2 


ewreeensw woe Pp fee wcoewooore 


Hua 
Dror 


od and awd ed 


R2=R2+DEL*RG 
ADsA2+DEL¥AG 
T2=22+DEL*TU 


(act) =asp*ra | 
~T22A24P1 
Rf) =R 1-R2 
“XM (I) =KM(I-1) ¢DEL#HSP*cos (T2) | 
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=XT #DEL*TSP*cos 
=YT +DEL*TSP*sin 
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MTA= MAL#LOS 


[ae | 


=0 


Sued £402) 33 
E 7 {tr ants #sen (PI 1 | 
i 


90 


Re eal ee oer ee 


Call SCREEN 


matter 


XM (I) =XM (I-1) #8 SPB*cos (T2) *DEL 
YM (I) =YM (I-11) +MSP¥*Sin (T2) * DEL 
XT (I) =XT (1-1) +TSP*cos(T1) *DEL 
YT (I) =YT(I-1) +TSP*sin (T1) *DEL 
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DO=NA V4 SP*TAC#COs (41) 78 (1-1) ------ 
D9=A (I-1) * (2*R2+NAV *MSP¥ cos (AZ) ) /R (I-11) 
[ a3=08-p9 | 


“D2=(TSP*sin (Al) -"SP¥sin (A2)) /R(I-1) 
[ (PSPesin (orac/TSB RG) | 


D7=NAV*D2 
P2=D6-D2 
_P3=07-D2_| 
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‘ R(I) =R (I-1) +DEL*R2 

A(T) =A (IL-1) #DEL*A3 

: P1=81+DEL*D2 

7 T1=f1+DEL*D06 

| T2=22+DEL*D7 
A1=A 1#+DEL*P2 

A2=A2+DEL*®P3 
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T5=DEL* (I-1) 
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Call SORT 


[ Call OUTPUT ce) 
J “Call SCREEN 


Do you want to run the 
same problem again? 
}ves 
LOS=RATODE (LOS 
TAL=RATODE(TAL 


MAL=RATODE (MAL 
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APPENDIX B. TRAJECTORY MODELS PROGRAM LISTING 


Following the next page is the program listing for the 
Trajectory Models program. [t consists of three main divi- 
sions; the executive routines, the FORTRAN IV computational 
program, and the FORTRAN IV plotting program. The 2xecutive 
routines are used to establish required file definitions, 
initiate operation of the somputational program, supervise 
the transfer of data to the plotting program, and provide 
operational information to the program user at appropriate 
times. 

The computational program, LPATH FORTRAN, consists of 
four subprogram divisions. The MAIN program accepts the 
input data and performs th2 calculations for the line of 
sight and the proportional navigation problems. Subroutine 
SCREEN is used to prompt the user to clear the terminal 
screen for proper positioning of the output. Subroutine 
SORT determines the largest missile acceleration value and 
the value ranges of the X and Y position coordinates for 
plotting reference. Finally, subroutine OUTPUT formats the 
calculated data and provides it to the user, the printer 
file, and the plot program data file. 

The plot program, PATHPLOI FORTRAN, is structured for 
direct submission to the MVS batch reader from the 3278 ter- 
minal. No carde need to be punched or read. Attention must 
be given to the correct JCL accounting data in the first 
line of the program; those shown in the listing are for 
illustrative purposes only. fhis program receives the data 
from the computational program and produces a singl2 geo- 
graphic chart of the encounter in the encounter ovlane. Mul- 
tiple problemas, up to ning will overlay on the single 
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APPENDIX C. WARHEAD DESIGN PROGRAM FLOWCHARTI 


XDEN,GO,VD,CDEN, RKILL 


Gc 
D,LT, 47, PKH ALT ,TTH, MA 


Present input data to user 


oe 2ee we oe eee ow p Men eo een ee ww oe 


>C=GC*GC 
T=LT*12 
WI=WT*12 


ALI>80000? 


Sa cea aE eS dels ESS ining CREE OER SI alates. ree ae Die Veni  - e  CPie ath ar 


I TE i at a a a, “Ler 
Py if é 
chee 3 


y 
t 


[ 


LWH=LTD*DW 
BETAS Seat lie 2 yal See 
RCRIT= HR) 7 (2¥*t 


(LT-Li 


RKILL>RCRIT? 


IFRC R=0 
A=PKH*WI/ (2 *PI*®RKILL) 


Ey ead tone ene EOBSEal oc ek ee eee he 
A=PKH* WT*LT/ (2* PI*RRILL* ( LWA 2*RKILL*tan (BETA) )) 
~CEM= (VI*VI/GC) /(1-(VI*VL/(2*GC) ) ) 


Write headier for PD table 


“gies aitee TTT] 


[~cali sizer 


oe ew @ 2 2e Bee Oe Be eee ee Re@eew west ae ewe nweaooewa oe 


Print PD(I),WHW,CTH, DI, FRAGS, HITS, PRLEN 


IFRBIG=1? 


“frite message indicating 
maximum value for Pi 


109 


wee dee em 


Naa ma 


eau. 


a4 ee 


i ee 


. 
7 
Sk ee A es ee 


Seon? We 


[ Input PDF ] 


=(VI*VI/GC) /(1- {VI *V 
CTM=(VI*VI/ ab ({) 283 1/(2*GC))) 


“Output WHW, XW,2W¥,CTH;DI,FRMLB(1),_ 
[ PREBN CP RAGS WETS. PDI) OC iy CRORES | 
“Gc=GCe*, 5 
LI=LT/12 
WI=WT 712 


pcobl em? 


[ana] 


110 


eae re en 


ee ee ee eS ‘hale ane a ha ¥ nab athe t0e a ee yaad 


es” 


R= (TTH*MAL) 7 (PRAREA**, 


cee SB i Salk oes AE ot 


[ PRPRASSEROER STDGGBER 227) | 


BLY (J) +1,.925 46x (1431,68754R) - (564. 1957*R*R) 
+(136. 706 sper a) (837 4O7 Rey 
BLV (J) =BLYV (J) *(. 26/CDEN) 


K=RHO*®*PRAREA*0.65/PRMLB(J 
V(J) =VHIT(I eoanan L 
M=V¥ (J) /SOOND 


‘ T=(1+.2*M**2) *TAIR 


Snuibis Haniiiiiaana 


s [ Return ] 


22 wee Bee ee eee e eee 


[ Subroutine SIZER 


IFRBIG=0 


NO 
BU= (2*CTN* (CDEN/EDEN) + 1) #FRALB( 1) *RKILL*12 
= 58C DEN *OKH*WT 
RW= ((-BU/BL) *alog (1-PD (I) )) ** (1/3) 


ee 2 ewe @e wee ee Oe ewaee fog GB e ee woe eseoeoees 2 ee weeoewe 


RW<DW/2? 


2 eee Sela 99. 
ee 


lll 


none hs SAE ROR T ny, I AAT SI te eS TSE tae NE Nite Na a tld Aes Ania Mihi eee 


shaal aca bot ba 
BL=5*C 


aati SS oN 


el ae ea 


~TPRBIG=1____ 
PD (I) =PD (I) *.99 
DI=0.0 


Tape tage 
a te 


GAMMA=1~ (1-1/8) **.5 


Bren ty lea 
CTH=GAMMA*R 


GK ee he 


é 


"_ xWe (PI® (RW-2TH) #*2) *LWH*X DEN 
CW=PI*(R Wa RW- ees © TH) **2) #LWH*CDEN 


CTH=R aW- (RWARW ~XW/ (PI *#LYH*CDEN 
2-CTH) ##2 XW / (XD EN*PI 


32 DI /DW 


CTN= WTinases a (tf eee day | 


CIN-C TM] <.00001TP >> 


cCTM=CTN 


112 


IN or Ake oe eae el ch Oh EG 


FRAGS= Sos bet hl A 
LEN= (FR 


on nn 


wae . 
—' 


See es es 


. *. ems Sap ie 7 ‘alee 5 feet x fi : : x 1 Y . es - é - 
a Re ans = : stele on meg IT aE NP gk CEN Sate TE VO Ie 2 e865 ei hal lle 256 Mbe Sa 


APPENDIX D. WARHEAD DESIGN PROGRAM LISTING 


This program has two major sections; the executive rout- 
ine and the PORTRAN IV comoutational program. The executive 
routine establishes the reuired file definitions and initi- 
ates operation of the computational program. 

The computational program, LBOMB FORTRAN, consists of 
four subprogram divisions. [he MAIN program accepts the 
input data, calculates the atmospheric characteristics, and 
formats and displays the output to the user and sends it to 
the printer file. Subroutine VEL calculates the initial 
velocity required to propel a given mass a specified dis- 
tance through the atmosphere with a particular residual vel-~ 
ocity remaining. It also iatermines the ballistic limit 
velocities for the situation. Subroutine SIZER sizes the 
warhead for a given Pd value. It also produces the charge- 
to-mass ratio, the number of fragments, the fragment size, 
and the average number of hits received by the target. 
Subroutine SCREEN prompts the user to clear the terminal 
screen for proper positioning of the displayed data. 
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APPENDIX F. PROPULSION SIZINS PROGRAM LISTING 


Following this page is the program listing for the Pro- 
pulsion Sizing Program. It has two segments; the 2xecutive 
routine and the FORTRAN IV computational program. The exe- 
cutive routine establishes the required file definitions and 
initiates operation of the computational program. 

The computational program, LPROP FORTRAN, consists of 
six subprogram divisions. The MAIN program accepts the 
input data from the user and performs the guiding calcula- 
tions for the booster and sustainer motors. Subroutine 
PRATIO determines the ambient pressures at the design alti-~ 
tudes. Subroutine PCRAT defines the optimum chamber: pres- 
sure to ambient pressure ratio with respect to the case 
material properties. Subroutine ARAT solves for the area 
ratio of the nozzle and tries to size the nozzle to the mis- 
sile diameter by varying the chamber pressure, characteris- 
tic velocity, and thrust coefficient. Subroutine SCREEN is 
used to prompt the user to clear the terminal screen for 
proper positioning of the output. And subroutine OUrPOT 
formats the calculated solutions and provides them to the 
user and to the printer file, if so directed by he user. 
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APPENDIX G. AERIJDYNAMIC COEFFICIENTS PROGRAM LISTING 
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This program is divided into three major subdivisions; 
the executive routines, the FORTRAN II/IV computational pro- 
gram, and the FORTRAN IV plotting program. The executive 


routines are used to establish the required file defini- 
tions, initiate operation of the computational program, 
supervise the transfer of jiata to the plotting program, and 


aire x 


PEER the, arscncenictiee 


provide operational information to the user at appropriate 


times. 
The computational program, LAERO1 FORTRAN, consists of 


oq. aoe 
PY Sorcerer ees og 


five subprogram divisions. [he 4AIN program accepts the 
input data, conducts calc lations not done by other subrout- 
ines, formats the calculated data and provides the output 
:! data to the user, the printer file, and the plot program 
data file. Subroutine GEOS0B calculates the missile wetted 
area and the Reynolds numbei per foot based on the flight 
altitude. Subroutine CLASJB calculates the aerodynamic sur- 
face lift-curve slopes. Subroutine CATSUB calculates center 
of pressure locations, cross-flow drag coefficients, and 
interference factors. Subroutine SCREEN prompts the user to 


awe 


: clear the terminal screen for proper positioning of the out- 
| put. 
{ The plot program, AEROPLOI FORTRAN, is structured for 
direct submission to the MVS batch reader from the 3278 ter- 
| minal. No cards need to ba panched or read. However, 
attention must be given to the correct JCL accounting data 
| in the first line of the program; those shown in ths list- 
ing are for illustrative purposes only. This program 
receives the data from the computational program and pro- 
duces a group of six charts for each Mach number entered. 
ee The plots represent the relationships 2£ Cl to alpha, Cm, 
cn, Cd, Ca, and Cdi. The program can produce up t) 24 sets 
of plots for a single run. 
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FILE: LAEROI EXEC a NAVAL POSTGRADUATE SCHOOL 
FILEDEF 08 DISK LAERO? OUPPOL AD (RECFM VA BLOCK 133 PERM 
ee DISK LAERO1 PLIT AO (RECFM VA BLOCK 80 PERM 

YOU WILL HAVE THE OPPORTUNITY TO OBIAIN BOTH A HARDCOPY 
PRINTOUT AND A SET OF PLOTS POR ONE SET OF INPUT PARAMETERS 
EACH TIME YOU ENTER THIS PROGRAM. THE PROGRAM MAY BE 
RE~RUN CONTINUOUSLY AND YIO WILL HAVE THE OPTION f9 CHANGE 
INPUT PARAMETERS FOR EACH SUCCESSIVE RUN, BUT YOU CAN 
OBTAIN THE PRINTOUT AND PLOTS PERTAINING [TO THE LAST RUN 
ONLY. IF ADDITIONAL OOTPJT LS REQUIRED, RE-ENTER THE 
PROGRAM. 

&END 

LOAD LAERO1 

START 

& BEGTYPE 

TO OBTAIN THE HARDCOPY PRINIOUT OF [HE DATA TABLES, TYPE 
AND ENTER: 


LAEROIPR 
hed THE VERSATEC PLOT OF THE TABULAR DATA, TYPE AND 


LAEROMVPL 
SEND 


FILE: LAEROPR EXEC A NAVAL POSTGRADUATE SCHOOL 


eaeenee LAERO1 OOTPUT (LINEZSIUN 70 


PE 
THE OUTPUT WILL BE PRINTEO ON THE VM PRINTER IN ROOM 140 
AND WILL BE IDENTIFIED WITA YOUR USER NUMBER AND LAST NAME. 
tans REQUIRES SEVERAL MINUTES TO OBTAIN THE PRINTOUT. 
PILE: LAEROWWL EXEC A NAVAL POSTGRADUATE SCHOOL 
COPY LAERO1 PLOT A LAERO PLOTDATA A PLOT LAERO? A 
EKEC SUBMIT PLOT LAERO1 A 
ERASE PLOT LAERO? A 
&BEGTYPE 
THE BLOt WILL BE DRAWN IN THE COMPUTER ROOM AND PLACED OVER 
THE ALPHABETIZ OUTPOT FILE IN ROOM 140. IT WILL BE 
IDENTIFIED BY TABROPILE" AND USUALLY REQUIRES MANY MINUTES 
faeces DAYS) TO OBTAIN 
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